Ordering and Spin Waves in NaNi02 : A Stacked Quantum Ferromagnet 
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Neutron scattering measurements on powder NaNi02 reveal magnetic Bragg peaks and spin waves 
characteristic of strongly correlated s=l/2 magnetic moments arranged in ferromagnetic layers which 
are stacked antiferromagnetically. This structure lends itself to stacking sequence frustration in the 
presence of mixing between nickel and alkali metal sites, possibly providing a natural explanation 
for the enigmatic spin glass state of the isostructural compound, LiNi02 • 

PACS numbers: 75.25.+Z, 75.40.Gb, 75.40.-s 



I. INTRODUCTION 

Low dimensional quantum magnets are well appreci- 
ated as fertile ground for spin liquid and other exotic 
quantum states of matter^. NaNi02 and LiNi02 are 
isostructural nickel-based quantum magnets with layered 
triangular structures. The enigmatic magnetic phase be- 
haviour associated with LiNi02 has been the subject of 
speculation for two decadesSiiiSiSiLSiS. This speculation 
has been largely fueled by the absence of observed transi- 
tions to long range magnetic and orbital order in LiNi02 
. In contrast NaNi02 is known to undergo a Jahn- Teller 
structural distortion at high temperatures signalling or- 
bital order^j and a low temperature anomaly in its sus- 
ceptibility is often associated with long range antiferro- 
magnetic order—. However, like LiNi02 , no definitive 
magnetic neutron scattering signature of such an ordered 
state has been observed in NaNi02 until now. We report 
the observation of magnetic Bragg peaks and spin wave 
scattering in NaNi02 , which establishes its relatively 
simple magnetic structure and which may shed light on 
why LiNi02 has difficulty finding an ordered state at low 
temperatures. 

Both NaNi02 and LiNi02 are comprised of stackings 
of triangular layers which alternate between NiO and 
AO, where A is either Li+ or Na+, as shown schemati- 
cally in Fig. 1. The distance between neighboring trian- 
gular layers of magnetic Ni along the c-direction is about 
5.2 A. Most of the extensive preceding work associated 
with the magnetism in these materials, assumes a low 
spin Ni'^"'", s=l/2 magnetic moment. However some dis- 
cussion of nickel in the s=l Ni^+ oxidation state bound 
to an s=l/2 hole on a neighbouring 0~ ion, resulting 
in an s= 1/2 Zhang-Rice singlet has been offerediSii^. In 
both scenarios, either the nickel spin itself, or the com- 
posite NiO spin is an extreme s=l/2 quantum mechani- 
cal entity, and a saturation magnetization corresponding 
to 1 per formula unit is observed in both NaNi02 and 



LiNiOjii. 



The possibility of antiferromagnetic coupling between 
these quantum moments within the triangular plane has 
generated much interest in the materials. The potential 
for geometrical frustration and exotic quantum ground 
states under such conditions have been well appreciated 
since Andersoni^ suggested a collective singlet, resonat- 
ing valence bond ground state without long range order 
for such a system. Indeed this possibility figures promi- 
nantly in the discussion of both early and more recent 
experimental work on LiNiO^SsI,. More generally, it is 
well known that the combination of antiferromagnetism 
and certain lattice symmetries based on triangles and 
tetrahedra leads to phenomena known broadly as geo- 
metrical frustration—. 



As mentioned above, the possibility of such exotic be- 
haviour has been motivated, in part, by the absence 
of direct experimental signatures of conventional mag- 
netic ordering, such as the observation of magnetic Bragg 
peaks by neutron diffraction. Given the difficulties as- 
sociated with the study of these materials, it is small 
wonder that progress has been slow. To date they have 
been available only as poly crystalline materials, and with 
a single s=l/2 magnetic moment per formula unit, the 
moment density is low. In addition, LiNi02 contains Li, 
whose strong neutron absorbtion is problematic for neu- 
tron scattering measurements. Finally, it is known that 
the similarity in ionic radii between Li"*" and Ni'^+ , both 
^ 0.7 A, leads to mixing between these two sublattices'^', 
further complicating the magnetic behaviour of LiNi02 
. With an ionic radius of ~ 1 A, Na+ is a much larger 
ion, which makes such mixing unlikely. For these latter 
two reasons, one expects NaNi02 to be more amenable 
to a neutron scattering study than LiNi02 . 
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FIG. 1: (a) The schematic magnetic structure of NaNi02 be- 
low Tjv ~ 23 K. (b) Stacking sequence frustration induced by 
impurity Ni spins on the alkah metal sublattice for antiferro- 
magnetic (J" > 0) coupling, (c) The same stacking sequence 
is favored by ferromagnetic (J" < 0, bottom) coupling. 



II. MATERIALS AND EXPERIMENTAL 
METHODS 

NaNi02 crystallizes into the rhombohedral space 
group R3m at high temperatures, before undergoing a 
cooperative Jahn- Teller distortion leading to a structure 
with the space group C2/m below ~ 480 KiflJi. This 
structural phase transition lifts the orbital degeneracy 
between the |3z^-r^> and |x^-y^> states within the Bg 
doublet in NaNi02-''. The resulting room temperature 
structuroiS in NaNi02 is characterized by lattice param- 
eters a=5.31 A, b=2.84 A, c=5.57 A, and ^=110.44 de- 
grees, and the structure no longer displays edge-shared, 
equilatoral triangles within the a-b plane. 

A similar structural distortion does not obviously oc- 
cur in LiNi02 although anomolies in the susceptibility 
are reported near 480 K in LiNi02^, which may be re- 
lated to a local structural distortion that cannot attain 
long range order due to mixing of the Li and Ni sublat- 
tices. 
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FIG. 2: top panel: Field-cooled and zero field cooled suscep- 
tibilities of NaNi02 are shown. Bottom panel: The temper- 
ature dependence of the integrated intensity of the (0,0,1/2) 
magnetic Bragg peak at Q~ 0.6 A is shown, indicating 
a continuous transition near Tjv 23 K. The inset shows 
net scattering at two different temperatures below Tjv as de- 
scribed in the text 



The lack of obvious orbital ordering in LiNi02 and 
the orbital degeneracy which would arise in its absence 
has motivated many theoretical authors to consider cou- 
pled orbital and spin degrees of freedomSiSiSiiSiiSiSSiSi as 
a mechanism to suppress magnetic ordering at low tem- 
peratures. Other recent work suggests orbital and spin 
degrees of freedom are decoupledi^. 

Stoichiometric amounts of Na202 and NiO were mixed 
and pelletized in an Ar atmosphere. These pellets were 
subsequently annealed at 973 K in O2 for 70 hours with 
one intermediate grinding. We prepared a 30 gram poly- 
crystalline sample of NaNi02 in this manner and charac- 
terized a small part of it with SQUID magnetic suscep- 
tibility techniques. The resulting field-cooled (FC) and 
zero field cooled (ZFC) dc susceptibilities are shown in 
the top panel of Fig. 2. The magnetic phase transition 
near Tn^ 23 K is immediately clear as a peak in x(T). 
In addition, a break between the FC and ZFC suscepti- 
bilities below '--^ 10 K is seen, indicating some glassiness 
within the ordered state at these low temperatures. 

Time-of-fiight neutron scattering measurements were 
performed on this 30 gram sample of NaNi02 using the 
Disk Chopper Spectrometer (DCS) at the NIST Cen- 
ter for Neutron Research. The DCS uses choppers to 
create pulses of monochromatic neutrons whose energy 
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FIG. 3: Maps of neutron scattering from NaNi02 at three 
temperatures near and below Tjv ~ 23 K; at a) 5 K, b) 10 K, 
and c) 22 K. The bottom right panel, d, shows the results of 
a linear spin wave theory calculation, at 10 K which can be 
compared directly to the experiments. The linear intensity 
scale has been chosen to highlight the inelastic spin wave 
scattering. 



transfers on scattering are determined from their arrival 
times in the instrument's 913 detectors located at scat- 
tering angles from -30 to 140 degrees. Measurements 
were performed with 3.2 and 5.5 Aincident neutrons. 
Using 5.5 Aincident neutrons, the energy resolution was 
0.075 meV^. 



III. EXPERIMENTAL RESULTS 

Typical data sets at low temperature are shown in Fig. 
3a, b and c, with the intensity scale chosen to highlight 
detail in the inelastic scattering. Cuts through the elastic 
scattering (integrating in tiuj between ±0.1 meV) of this 
data are shown in the inset to the lower panel of Fig. 2. 
Here we clearly see the appearace of new Bragg peaks 
below Tn ~ 23 K. These particular data sets correspond 
to the subtraction of a high temperature, 30 K, data set 
from data sets at the indicated low temperatures and 

are centered around the Q=0.6 A region of reciprocal 
space, the position at which the lowest-Q Bragg peak 
appears. 

This and several other very weak magnetic Bragg 
peaks 10~^ of the strongest nuclear Bragg peaks) are 
evident and can be accounted for both in position and 
relative intensity by a simple magnetic structure in which 
the s=l/2 moments on the nickel sites are arranged in 
ferromagnetic sheets within the triangular planes, and 



FIG. 4: Cuts through the maps of scattering at 10 K shown 
in Fig. 3b. These data approximate constant-Q scans at 
Q=0.52 A and 0.86 A ^ and clearly show a largely dis- 
persionless excitation at fiui ~ 0.7 meV, in addition to the 
Goldstone mode. 



the sheets are antiferromagnetically stacked along c, such 
that moments on neighbouring planes are tt out of phase 
with each other. The details of the modeling of this elas- 
tic magnetic scattering are consistent with a rather large 
ordered moment, close to 1 /is/Ni, and do not support 
composite Zhang-Rice singlet magnetic moments associ- 
ated with NiO. 

The Bragg peak shown in the bottom panel of Fig. 2 
is indexed as (0,0,1/2) and arises due to a tt phase shift 
across the 5.2 Abetween triangular Ni planes, thereby 
appearing at Q=7r/5.2 A'--^ 0.6 A . The integrated in- 
tensity of this Bragg peak is shown in the bottom panel 
of Fig. 2. Both this intensity and that of other magnetic 
Bragg peaks tend to zero near T^v^ 23 K, consistent 
with a continuous phase transition. The rounded nature 
of the order parameter near T^r is unusual, and possibily 
due to dimensional cross over. 

The details of the magnetic structure, to be presented 
elsewhere, are determined by the relative intensities of 
the measured magnetic Bragg peaks at low tempera- 
tures. Here we simply note that the form of the magnetic 
neutron scattering cross section requires that Bragg in- 
tensity at (0,0,1/2) be due to components of ordered mo- 
ment within the triangular plane, and that the relative 
intensities require that the moments make a relatively 
large angle with respect to the triangular plane. 

The inelastic scattering in Fig. 3 clearly shows spin 
wave excitations going to zero energy at the (0,0,1/2) 
magnetic zone centre, and reaching a zone boundary en- 
ergy of ~ 0.7 meV at Q~0.3 A . In addition, rather 
weak but clearly observable inelastic scattering is seen 
in a dispersionless inelastic feature at ~ 0.7 meV which 
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IV. LINEAR SPIN WAVE THEORY AND 
COMPARISON TO EXPERIMENT 



We have also carried out linear spin wave theory cal- 
culations appropriate to s=l/2 moments interacting on 
the Ni^+ sublattice of NaNi02 , with the following mi- 
croscopic Hamiltonian: 



|Q|(A ') 



FIG. 5: Maps of calculated neutron scattering from NaNi02 
using linear spin wave theory is shown at 10 K, and for four 
different easy plane orientations relative to the triangular 
basal plane. a=0, shown in a), corresponds to a coincidence 
between the easy plane and the triangular plane, while a=90, 
in d), corresponds to these two planes being normal to each 
other. The calculation for q=40 degrees in panel b) is the 
same as that shown in Fig. 3d) . These calculations show that 
a relatively large angle a is required between the magnetic 
easy plane and the triangular basal plane, in order for the 
dispersionless excitation near 0.7 meV to be observable. 



extends across all Q measured. This is most clearly seen 
in cuts through this data, which approximate constant- 
Q scans, and which are shown for two different Q's in 
Fig. 4. All of these excitations can be seen on both the 
neutron energy loss (-|-ve side) and neutron energy gain 
(-ve side) of zero energy transfer, as required by detailed 
balance. The inelastic data in Fig. 3, has been corrected 
for detector efficiency only. 

The inelastic scattering meaurements near T^v show 
the spin wave band along c* to soften completely, such 
that the envelope of the spin wave dispersion fills in with 
inelastic intensity, as shown in Fig. 3c. These measure- 
ments indicate the nature of the transition at Tn ~ 23 K 
is a loss of registry between well-correlated ferromagnetic 
triangular planes. This likely also explains the previous 
observation of small angle neutron scattering in LiNi02^ 
as originating from the collapse of the inelastic spin wave 

scattering, which occurs below Q~0.6 A . Such inelas- 
tic scattering would be integrated up in a diffraction, 
SANS experiment producing a temperature-dependent 
signal at small- Q, whose intensity peaks near Tn. 



n ^ jY, SiSj + X ^ s!s^ + SiSj (1) 

where the exchange integral within the triangular ah 
basal plane is ferromagnetic, and relatively strong, J=- 
2.5 meV, and that along the stacking c direction is an- 
tifcrromagnetic and relatively weak, J'=-|-0.16 meV. In 
addition we have included an anisotopy term, K, whose 
strength is set equal to J', and whose purpose is to pref- 
erentially restrict the spins to a plane perpendicular to z. 
All interactions are nearest neighbour only. The inclu- 
sion of this anisotropy term is phenomenological, produc- 
ing two transverse branches in the spin wave spectrum 
along c*, as is observed experimentally. Note that the 
Dzyaloshinski-Moriya interaction, which requires a lack 
of inversion symmetry between the sites whose coupling 
it mediates, is not allow by the space group of NaNi02 . 
Hence anisotropic exchange is expected to be the leading 
order anisotropic interaction. 

Imx(Q, Ti^) for such a system was calculated within 
linear spin wave theory at T=10 K and this result was 
angularly averaged to produce theoretical expectations 
appropriate to an inelastic neutron scattering study car- 
ried out on a polycrystalline sample. This is shown in 
Fig. 3d and in all the panel of Fig. 5. 

This calculation employs a z-direction in Eq. 1, nor- 
mal to the easy plane, which defines the easy plane. In 
general, the resulting magnetic easy plane makes some 
angle, a, with respect to the triangular basal plane. 
Therefore the magnetic easy plane is not neccesarily 
coincident with the triangular plane. The calculation 
shown in Fig. 3d employs an angle of a=40 degrees. Fig. 
5 shows the results of this calculation for the expected 
powder-averaged inelastic neutron scattering using four 
different values of a: a) a=0 degrees (easy plane coinci- 
dent with triangular plane), b) a=40 degrees, c) a=60 
degrees, and d) q;=90 degrees (easy plane normal to the 
triangular plane). 

As can be seen, a non-zero a is required in order for 
the dispersionless mode at ~ 0.7 meV to have observable 
weight. If the magnetic easy plane was coincident with 
the triangular basal plane, as they are at a=0, fluctua- 
tions out of this plane would be along c*, and not ob- 
servable along this direction, by virtue of the polarization 
dependence of the neutron scattering cross section. In 
this a=0 case, only the Goldstone mode is predicted to 
be observed. 
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As can be seen comparing Fig. 3b and 3d, the quan- 
titative agreement between the calculation of the spin 
wave spectrum and the measurements at 10 K is very 
good, although some minor discrepancies are evident. 
The measured spin wave excitations are identified as cor- 
responding to those propagating along the stacking di- 
rection. Furthermore, the zone boundary spin wave en- 
ergy, at Q=0.3 is given by Azs = 6S^J{J''^ + J' K) 

which reduces to IS.zb = GV^SJ' for the case where 
J'=K. Reading Azb —0-7 meV off the data, we get 
J'=0.16 meV. This direct measurement of J' can be com- 
bined with the measured Ocw~3.5 meV obtained from 
fits to the high temperature susceptibilit y^'-- to produce 
an estimate for the stronger in-plane ferromagnetic ex- 
change. Using kBQcw = -^^^f^(3(J + J') + K) we 
obtain J=-2.5 meV. 



V. DISCUSSION AND CONCLUSIONS 

Our measurements of the spin structure and dynam- 
ics in NaNi02 directly reveal a rather simple magnetic 
structure corresponding to an antiferromagnetic stacking 
of ferromagnetic sheets. Such a structure is very sensi- 
tive to frustration of the stacking sequence, and there- 
fore frustration in the ability of the structure to attain 
true three dimesional long range order, due to mixing of 
the alkali metal and transition metal sublattices. This 
is precisely the nature of chemical disorder shown to be 
relevant in LiNi02 at the 1-3 % leve&S. 

This scenario is illustrated schematically in Fig. 1. 
This drawing shows 3 sets of interacting ferromagnetic 
layers of Ni moments. Figure 1(a) shows the simple an- 
tiferromagnetic stacking of ferromagnetic layers deter- 
mined from our neutron measurements. Figure 1 (b) and 
(c) show the consequences of impurity spins which cou- 
ple either antiferromagnetically (b) or ferromagnetically 
(c) to neighboring Ni layers. It is clear that misplaced 
magnetic Ni ions sitting in the alkali metal triangular 
planes provide an exchange pathway along the stack- 
ing direction which will frustrate the stacking sequence 
independent of whether these impurity spins couple fer- 
romagnetically or antiferromagnetically to the moments 



on Ni triangular planes above and below. This is due 
to the fact that two such pathways link neighbouring 
Ni triangular planes, thus the sign of this interaction is 
irrelevant. Furthermore, as these impurity nickel spins 
are a factor of two closer along the stacking direction to 
spins within the Ni triangular planes than those residing 
within the Ni triangular planes, the magnitude of this 
frustrating Ni - impurity Ni - Ni interaction is expected 
to be relatively strong. 

If, as has been argued^ii^, orbital degeneracy is not re- 
sponsible for it, a natural and simple explanation for the 
spin glass phase in LiNi02 ensues as a consequence of 
this structure. Samples of LiNi02 display weak mixing 
between the Ni and Li sublattices, providing a sample de- 
pendent tuning of the strength of the stacking sequence 
disorder. This leads to a glass transition which depends 
on the precise level of disorder in the sample and which 
will occur at temperatures on the order of that associ- 
ated with phase coherence between strongly correlated 
ferromagnetic layers, Tjv ~ 23 K. This scenario is qual- 
itatively similar to the ferrimagnetic clusters previously 
proposed22iSi to explain the spin glass phase in LiNi02 , 
with the exception that the stacking sequence frustration 
occurs independent of the sign of the coupling between 
impurity spins and the magnetic layers. 

To conclude, new elastic and inelastic neutron scat- 
tering measurements have directly determined the sim- 
ple antiferromagnetic structure of NaNi02 below T^r ~ 
23 K. These measurements allow a microscopic under- 
standing of the magnetically-ordered state in NaNi02 
and may provide the basis for a simple explanation of 
the phase behavior exhibited in LiNi02 . 
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